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, Abstract. In this note we document the existence of a finitely generated rational cone 

^"*^ ■ that is not covered by its unimodular Hilbert subcones, but satisfies the integral Caratheo- 

dory property. We explain the algorithms that decide these properties and describe our 
experimental approach that led to the discovery of the examples. 
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1. Introduction 



o 

u 

Let C C Mr be a finitely generated rational cone, i. e. the set of all linear combinations 

a\X\-\ Ya n x n of rational vectors x\ , . . . ,Xd with coefficients from R + . We can of course 

assume that £ Z d , i= 1, . . . ,n. In this note a cone is always supposed to be rational and 
finitely generated. Moreover, we will assume that C is pointed: if x, —x £ C, then x = 0. 
^ ■ Finally, it is tacitly understood that C has full dimension d. 

The monoid M(C) = Cfl Z d is finitely generated by Gordan's lemma (for example, see 
[|2l Section 2. A]). Since C is pointed, M is a positive monoid so that is the only invertible 
element in M(C) . 

It is not hard to see that M(C) has a unique minimal system of generators that we call 
its Hilbert basis, denoted by Hilb(M(C)) or simply Hilb(C). It consists of those elements 
z 7^ of M(C) that have no decomposition z = x + y in M(C) with y, z ^ 0. 

We want to discuss combinatorial conditions on Hilb(C) expressing that C or M(C) 
is covered by certain "simple" subcones or submonoids, respectively. To this end we 
define a u-subcone of C to be a subcone generated by vectors x\,. . . ^xj £ Hilb(C) that 
form a basis of the group Z d . In particular, xi are linearly independent, and if 
just this weaker condition is satisfied, then the cone S generated by x\, . .. ,Xd is called an 
f -subcone. In this case we let T(S) denote the subgroup of Z d generated by x\ , . . . , x c \ and 
E(S) the submonoid of 7L d generated by x\, . . . ,Xd- Note that 5 is a w-subcone if and only 
if T(S) = Z d , or, equivalently, E(S) = S n Z d . 

One says that C satisfies (UHC) if C is the union of its w-subcones. The letter U stands 
for unimodular, H reminds us of the condition that the generators of the a-subcones be- 
long to Hilb(C), and C simply stands for cover. 

A weaker condition than (UHC) is the integral Caratheodory property (ICP). One says 
that C has (ICP) if every element of M(C) can be written as a linear combination of at 
most d elements xi £ Hilb(C) with integral nonnegative coefficients a,. The terminology 
is motivated by Caratheodory 's theorem: let y\ , . . . ,y m be a minimal system of generators 

of the cone C; then every element y £ C is a linear combination y = a\yi x H h a^yi d 

with nonnegative real coefficients. 

Both (UHC) and (ICP) can be formulated more generally for positive affine monoids 
McZ . However, it is easy to see that every monoid M satisfying (UHC) is given in 
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the form M = M + M n T, d . By a theorem of Bruns and Gubeladze Q3 Theorem 6.1] the 
same holds true if M satisfies (ICP), provided the group gp(M) generated by M equals 
7L d . In loc. cit. it is also shown that (ICP) is equivalent to the formally stronger condition 
that M is the union of its submonoids £(5). (This condition is called (FHC) in 0]|.) The 
equivalence is crucial for our note, and therefore we reproduce the statement and its proof 
in Theorem [2l 

While we view (UHC) and (ICP) as structural properties of (normal) affine monoids, 
these properties have first been discussed in the context of integer programming: see 
Cook, Fonlupt and Schrijver [[6| and Sebo flU. 

It was asked by Sebo [8] whether every cone C has (ICP) or (UHC), and he proved that 
(UHC) holds if d < 3. He actually proved a stronger statement: C has a triangulation by 
M-subcones. A counterexample to (UHC) in dimension 6, called C\q in the following, was 
found by Bruns and Gubeladze [1], and then verified to violate (ICP), too, in cooperation 
with Henk, Martin, and Weismantel [3]. Despite the existence of the counterexample, one 
can fairly say, at least heuristically, that almost all cones satisfy (UHC). 

It remained an open problem whether (UHC) is strictly stronger than (ICP). In this note 
we want to document the existence of cones that satisfy (ICP) but fail (UHC), explain the 
algorithms that decide (UHC) and (ICP), and describe the experimental approach that led 
to the discovery of the examples. 

All our experiments seem to indicate that C\q is the core counterexample to (ICP) and 
(UHC). In fact, all counterexamples to these properties that we have been found contain 
it. It would be very desirable indeed to clarify the situation in dimensions 4 and 5. 

Acknowledgement. The author is very grateful to Joseph Gubeladze for inspiring dis- 
cussions and to the Mathematisches Forschungsinstitut Oberwolfach where the first steps 
of this project were taken during a joint visit within the MFO's RiP program. 

2. Deciding UHC 

Let us say that x E C is u-covered if it is contained in a w-subcone. A subset of C is 
w-covered if each of its elements is w-covered. Using this simple terminology, we can 
describe an algorithm deciding (UHC); see Tabled] 

In the algorithm we use a function named split. It decomposes D along a support 
hyperplane H of a M-subcone U such that D fl H y ^ as well as D fl H < ^ 0. Such a 
hyperplane does indeed exist if int(D) flint(£/) 7^ 0, but D U. The cones produced are 
D\ — D n H + and D2 = D DH~ . (The open halfspaces determined by H are denoted by 
// > and H < , and H + and H are the corresponding closed halfspaces.) 

It is easy to see that the algorithm terminates: there are only finitely many hyperplanes 
by which we split subcones. Therefore only finitely many subcones can be created. 

In order to check the correctness of the algorithm, observe that at each step in the for 
loop in unicover none of the w-subcones U\, . . ., t/ ; _i of C intersects the interior of D. 
In fact, the index i is only increased if int(D) fl int(t/,) = 0. (In the recursive call the index 
i that has been reached at the parent level starts the for loop at the child level.) 

Thus, when the loop terminates with i = N+l, none of the w-subcones of C intersects 
the interior of D. So D, and consequently C, indeed contains a vector that is not u- 
covered. Conversely, if C contains such a vector x, then on each level of the recursion tree 
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UNICOVER(D,n) 

1 for i^ntoN 

2 do 

3 if DcU t 

4 then return 

5 ifint(D)nint(E//)^0 

6 then (D U D 2 ) «- sPLn(D,Ui} 

7 unicover(Di,/) 

8 unicqver(Z>2,/) 

9 return 

10 output ( D not w-covered ) 

1 1 return 

MAIN0 

1 Create the list U\ , . . . , U?j of w-subcones of C 

2 unicqver(C, 1) 

Table 1 . An algorithm deciding UHC 




one finds a subcone D containing x, and for such D the condition D C Ui can never be 
satisfied. Therefore there exists an end node of the recursion tree at which the loop is left 
with i = N+l. 

The pseudocode in Table[T]is a somewhat simplistic sketch of the actual implementation 
since it is necessary to cope with substantial memory requirements. For example, the list 
Ui, ...,Un is not produced a priori, but extended whenever necessary and always kept 
as small as possible. Moreover, all allocated memory is recycled carefully within the 
program. 

Instead of starting the covering algorithm with the full cone C, the actual implemen- 
tation uses the output of a preprocessor that computes several triangulations A\,.. .,A t 
of C, The input to uni cover (and also to caradec below) is the list of intersections 
D\ fl • • • HD t where Di is a nonunimodular simplicial cone in A;. 
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3. Deciding ICP 

Let us first fix some terminology that parallels that for (UHC). An element x E Cfl 1, d is 
f -covered if it belongs to one of the monoids E(5) where S is an /-subcone, and a subset 
of C is /-covered if each of its elements is /-covered. 

The lemma contains the basic criterion by which we can check that C is /-covered. In 
the theorem following it, we will then see that this property is equivalent to (ICP). 

Lemma 1. 

(a) Let G\ , . . . , G n be subgroups of 7L d , and let N be a residue class of 7/ modulo 
Gi n • • • n G n . Then N C Gi U • • • U G„ if and only ifNn {G\ U • • • U G n ) ^ 0. 

(b) Let Si,...,S n be f-subcones ofC, each containing the d-dimensional subcone D 
of C. If every residue class oflr modulo T{S\) fl - • • fl T(S„) meets T(5i) U • • ■ U 
T(S n ), then D is f -covered. 

(c) Let D be a d-dimensional subcone of C with the following property: for every 
f -subcone S either D C S or int(D) fl int(S) = 0. Furthermore let Gu be the 
intersection of the groups T(S), S D D, and Ho their union. Then D is f -covered 
if and only if every residue class of If 1 modulo Go meets Ho- 

Proof, (a) Suppose that N fl (G\ U • • • U G n ) ^ 0, and let x be an element in the intersection, 
x E Gi. The subgroup G' = G\ fl • • • fl G n is contained in G,-, and so N = x + G' C G,. The 
converse implication is trivial. 

(b) Let x E D fl II 1 . It follows from (a) that x E r(5,) for some i. But x E Si, too. 
Therefore x E T(S) flS,- = (At this point we use that the generators of 5,- are linearly 
independent.) 

(c) It only remains to show the necessity of the condition. For it we only need to 
observe that every residue class Af of Z d modulo Go meets int(D). By hypothesis on 
D, an element x E iV flint (D) is /-covered if and only if x E E(5) for some /-subcone S 
containing D. □ 

We include the next theorem and its proof for the convenience of the reader. It is 
a simplified version of ffi Theorem 6.1] whose proof contains the crucial ideas for the 
algorithm deciding (ICP). 

Theorem 2. Let M C Z d be a positive affine monoid such that gp(M) = If 1 . IfM satisfies 
(ICP), then M = M+M fl and every element ofM is f -covered. 

Proof. We dissect C = IR + M along all the support hyperplanes of the cones spanned by 
linearly independent vectors x\,...,Xd of Hilb(M) into elementary subcones. Set M = 
M+M f Z d and choose xEM. Suppose that x has no representation as a linear combination 

x = a\y\ H h ajyd, a\,...,ad E Z + and yi, . .. ,yd E Hilb(M) linearly independent. 

The element x belongs to one of the elementary subcones D, and as in the proof of the 
lemma it follows that there exists a finite index subgroup G of Z rf such that no element 

z of (x+ G) flint(D) has a representation a\y\ H h a^yd with a\,...,ad E Z + and 

yi, • • • ,yd £ Hilb(M) linearly independent. 

The crucial point is that M \M is contained in the union of finitely many hyperplanes 
(see [2, Section 2.B]), and the same applies to all elements of M that are linear combina- 
tions of linearly dependent elements of Hilb(M). But (x + G) fl int(D) is not contained 
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in the union of finitely many hyperplanes, and so must contain elements of M. This is 
impossible if M satisfies (ICP). □ 

caradec(Z), G, M, n) 

1 for i^ntoN 

2 do 



3 ifD c Si 

4 then M' <- 

5 for (xef, G/(Gnr(5 f ))) 

6 do 

7 ifx + y^r(5,-) 

8 then^' =^'u{jc+y} 

9 G<— Gnr(Sj), 

10 if^ = 

1 1 then return 

12 ifD £ Si and int(D) n int(S f ) ^ 

13 then (D\,D 2 ) <- split(D,5 ; ) 

14 CARADEC (Di ,G,M, i) 

15 CARADEC(D2,G,^,/) 

16 return 



17 output ( Z) not /-covered ) 

1 8 return 

MAIN0 

1 Create the list S\ , . . . , of /-subcones of C 

2 CARADEc(C,Z d ,{0},l) 

Table 2. An algorithm deciding ICP 



For the algorithm deciding (ICP) we have to enrich our data structure by those com- 
ponents that have shown up in the proof of the lemma. Subcones are replaced by triples 
(D, G,i?) where D is a subcone of C, G is a finite index subgroup of Z d , and S% is a list 
of residue classes in Z / G. In & each residue class is represented by a single vector that 
belongs to it, and in the algorithm (see Table |2]) the loop 

for (jcef, jeG/(Gnr(5;))) 

runs over all elements of ffl x G/(Gnr(5 ; )). 

Again it is clear that the algorithm terminates after finitely many steps: the number of 
hyperplanes that we can use to split subcones of C is still finite (though larger than for 
(UHQ). 

The crucial point for caradec is that at each step in the loop the /-cones S\,...,Sn 
satisfy the following conditions: 

(1) for each j<i- \ either Sj D D or int(S/) n int(D) = 0; 

(2) G is the intersection of all groups T(Sj), j < i — 1, for which D C Sf, 
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(3) is the list of those residue classes in Z rf /G that are not contained in the union 
of the groups T(Sj), j < i — 1, D C Sj. 
We have only to check that these conditions remain satisfied when 5 ; is tested against 
D. To this end let S^ , • • • , be those among Si , . . . , that contain D. 
If int(5;) PI int(D) = 0, then D (£_ Sj, and this case is done. 

If int(S;) flint(D) ^ 0, but D <£_ S\, then i is not increased (!) and all three conditions 
are inherited by both D\ and D2 '. among the Sj, j < i — 1, exactly Sk t , . . . , Sfc,„ contain D\ 
or D2, simply because Sj D £>i or Sj D D\ implies int(Sj) fl int(D) 7^ 0, and so Sj D D. 

But if Si D D, the bookkeeping is also correct. Evidently G is replaced by the correct 
group Gfir(S ; ). Next observe that all residue classes of Gnr(S ; ) that are contained in 
residue classes modulo G not appearing in M remain in T(Si Cl ) U • • • U r(5yt m ) . On the other 
hand, those that refine elements of must belong to r(S ; ) to be in T(Sk l ) U • ■ • U T(Sk m ) U 
r(5 ; ). The correctness of the algorithm follows now immediately from Lemma[TJ 

The biggest hurdle for it are the lists & of residue classes that usually become extremely 
long already in dimension 6. Moreover, along each branch of the recursion tree, several 
of them must be kept in memory. (This problem cannot be eliminated by a nonrecursive 
implementation.) 

The growth of the list & can be estimated. Set 

e = #(G/(G nr(5,-))) and e' = #(Z d /T(Si)). 

Each element x e & is involved in e vectors x+y. At most one of them lies in r(5,) since 
the vectors y belong to pairwise different residue classes modulo r(5,). Therefore 

#(^') > {e-\)#{M). 

If the elements of & are randomly distributed over the residue classes of Z rf modulo T(5 ; ), 
then the expected share of vectors x + y E T(Sj) drops to l/e'. 

Instead of keeping the lists of residue classes in memory, one could alternatively try 
to follow the recursion tree along the whole list 5i , . . . , Sn, compute only G along each 
branch and test the residue classes one by one only at the end nodes. However, this 
approach seems unfeasible since it derives no advantage from the case e = 1 , which for- 
tunately happens frequently and often stops the recursion before the end of Si , . . . , Sn is 
reached. 

The list 5i,...,5iv is actually scanned in growing order of the determinants of the Sj. 
This has turned out very effective, at least for those cones that satisfy (ICP). In fact, all 
cones in Table [5] with (ICP) are covered by /-subcones of determinant < 2. 

In addition to caradec we use a Monte Carlo approach for disproving (ICP). It reads 
the output of un 1 cover, computes a large number of vectors in the non-w-covered sub- 
cones of C and tests whether they are /-covered. 

Remark 3. caradec provides us with a precise measure for the failure of (ICP), namely 
the ratios #(^) /#(G) at the end nodes of the recursion tree. For the cone C\o (Table [3]) 
there is precisely one end node with M 7^ 0, and the ratio is 32/15552 = 1/486. The num- 
ber of non-/-covered vectors in the Monte Carlo test confirms the ratio rather precisely. 
For the cone C'( 5 (Table [5]) there is again a single group with & 7^ and ratio 

2, 4468, 480/2, 286, 144, 000, 000 w 1 .07/ 10 6 . 
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So the Monte Carlo test cannot be expected to be conclusive with < 10 6 test vectors. 

4. The search 

Let us recapitulate an important notion from [1J. An element x of Hilb(C) is called 
destructive if H' = Hilb(C) \ {x} is not the Hilbert basis of H + H'. We say that C is tight 
if every element of Hilb(C) is destructive. The crucial role of tight cones for (UHC) and 
(ICP) is illuminated by the following lemma [1 , Corollary 2.3]. 

Lemma 4. Let C be a cone that is a counterexample to (UHC) or (ICP). Suppose C is 
minimal first with respect to dimension and second with respect to #Hilb(C). Then C is 
tight. 

Remark 5. Updating the information in [1] we mention that tight cones exist in all dimen- 
sions d > 3. The first 3 -dimensional tight cone was found by P. Dueck. The smallest such 
cone found by the author has a Hilbert basis of 19 elements. The elements of the Hilbert 
basis in the extreme rays form a regular hexagon (with respect to the action of GL3(Z)) 
so that the cone has the dihedral group D(, as its automorphism group. The regularity is in 
indication that it may be the smallest possible tight cone. (Here and in the following the 
automorphism group of a cone C is always understood to be the automorphism group of 
the monoid CnZ d .) 

Our search for counterexamples has been based on the crucial Lemma @] We produce 
a set of random vectors, consider them as the generating set of a cone C, and then use a 
program named shrink to remove nondestructive elements of Hilb(C) until a tight cone 
is reached, (shrink is based on the same algorithm as normaliz; see BUO.) Almost 
always, C shrinks to the 0-cone, but sometimes a nontrivial tight cone emerges. Then 
un i cover, and possibly caradec, are invoked. 

When we started the search in spring 1998, we used cones over randomly generated 
lattice parallelepipeds. In May 1998 the search stopped with the counterexample Ciq. Its 
Hilbert basis is shown in Table |3] The cone Ciq has 27 support hyperplanes. 



Z\ = 


(0,1,0, 0, 0, 0) 


Z6 = 


(1,0,2,1,1,2) 


Z2 = 


(0, 0,1,0, 0, 0) 


Zl = 


(1,2,0,2,1,1) 


Z3 = 


(0, 0, 0,1,0, 0) 


zs = 


(1,1,2,0, 2,1) 


Z4 = 


(0, 0, 0, 0,1,0) 


Z9 = 


(1,1,1,2,0, 2) 


Z5 = 


(0, 0, 0, 0, 0,1) 


ZlO = 


(1,2,1,1,2,0) 



Table 3. Hilb(Cio) 



The reader should note that for the questions considered in this note we can always 
replace a given cone C by 0(C) where is an arbitrary transformation in GL^(Z). In this 
sense, C stands for a class of cones that are isomorphic under an integral isomorphism of 
M. d . We express this fact by speaking of different embeddings of a cone C. 

While un i cover showed that C\q fails (UHC), it was then verified in cooperation with 
Henk, Martin, and Weismantel that C\q also fails (ICP). (caradec was not written before 
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September 2006.) See Bruns and Gubeladze [1] and Bruns et al. [3] for more information 
on C w . 

The automorphism group of C\q is remarkably large: it is the Frobenius group F20 of 
order 20, which acts transitively on zi, ■ ■ ■ ,zio- (7*20 is the semidirect product of Z5 with 
its automorphism group Zg = Z4.) From the embedding above one can see that at least the 
dihedral group D5 C F20 is acting on C\q. All the remaining 10 automorphisms have order 
4 and swap {zi, . . . ,15} with {z 6 , . . . ,zio}. Moreover, z\,...,zw all lie in the hyperplane 

given by —5^1 + H h ?6 = 1- The convex hulls of {zi,...,zs} and {z 6 • • • ,Zio} are 

both simplices of dimension 4. 

Remark 6. It was communicated to us by F Santos that the lattice polytope spanned 
by Hilb(Cio) is a projection of the Ohsugi-Hibi polytope [7J. The projection leads to 
the following description of C\q fl Z 6 . Consider the complete graph £5 and decompose 
it into 2 cycles of length 5 as shown in Figure [2l Now choose the incidence vectors 




Figure 2. Cycle decomposition of £5 

(1, 1,0,0,0) etc. of the edges in the first cycle and prefix them with 0. Then prefix the 
incidence vectors (1,0, 1,0,0) etc. of the second cycle with 1. The resulting 10 vectors in 
Z 6 generate a monoid M isomorphic with C10 fl Z 6 . While this description is even more 
aesthetic than the one in Tabled it has the disadvantage that gp(M) is of index 2 in Z 6 . 

In the summer of 1998 a second counterexample C\i to (UHC) and (ICP) emerged. It 
has a Hilbert basis of 12 elements. We continued the search for two more years. The frus- 
trating outcome was that C10 appeared over and over again, but no new counterexample 
showed up (and even C12 did not return until November 22, 2006). 

The project was taken up again at the end of 2004 when our department had installed 
a dual processor Opteron system with very fast integer arithmetic. Nevertheless, the out- 
come of the search remained as disappointing as it had been before. 

Finally, in August 2006 we did what should have been done long before, namely com- 
pare C\2 with C\q\ it turned out that Hilb(Ci2) (in an embedding that had to be found!) 
extends Hilb(Cio) by two vectors. Relative to C\q, this finding explained why C12 fails 
(UHC) and (ICP), too: the extra a-subcones and /-subcones are not sufficient to cover all 
integral vectors in C\$. (It also shows that one cannot speed up shrinking by removing 
two vectors at a time.) 

However, this not very surprising a posteriori insight made it suddenly clear that there 
might be many interesting objects in the vicinity of C\q. Especially, when we approach 
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Cio along a shrink path, why should the stronger property (UHC) not be lost before (ICP)? 
After a modification of shrink we also applied uni cover to the, say, 6 last non-tight 
approximations to Cio, and within hours many new non-(UHC) cones emerged. Several 
of them defeated all Monte Carlo attacks on (ICP). It became clear that caradec had to 
be implemented, and it indeed recognized many non-(UHC), but (ICP) cones. 

Ironically, within a few weeks after we had given up our narrow-minded insistence 
on checking only tight cones, two new non-(UHC) such cones surfaced, both of them 
satisfying (ICP). They appear as C' u and C15 in Tabled Since all these cones contain 



Ci 2 : z' n = (2, 2, 1,4,1,3) C 15 : m = (2, 1, 0, 5, 1, 5) 

4 = (2, 3, 1,4, 1,2) w 2 = (1,0, -1,4, 0,4) 

w 3 = (0, 0,-1,1,0,1) 
C' n : z'l x = (0,-1, 2, -1,-1, 2) w 4 = (2,l,2,3,2,4) 

z'/ 2 = (l,0,3,0,0,3) w 5 = (1,1, 0,3, 1,2). 

Table 4. Additional vectors in Hilb(Ci 2 ), Hilb(C / 12 ), Hilb(Ci 5 ) 



Cio, we list only the extra vectors that complement the Hilbert basis of Cio (using the 
embedding given in Table [3]). The numbers of support hyperplanes are 39 for C12, 40 for 
C[ 2 , and 36 for C15. 

The most interesting after Cio undoubtedly is C' 12 , not only because it satisfies (ICP). 
Its automorphism group - certainly invisible from the embedding given - is again the 
Frobenius group F20. It is clear that F20 cannot act transitively on Hilb(C' 12 ), which rather 
decomposes into an orbit of 10 elements and one of 2. However, this is by no means an 
extension of the action of F20 on Cio since the orbit of two elements is {21,25}! Only 
a subgroup isomorphic to Z 2 x Z 2 restricts to Cio, showing that there are 5 conjugate 
embeddings of Cio into C' 12 , and each of them contains 21,25. 

The convex hulls of {21, . . . ,25,2^} and {26 ... ,210,2^} are both bipyramids over a 
tetrahedron. The bipyramids are situated in the parallel hyperplanes with the equations 
£i=0 and £i = l. Both C12 and C[ 2 have their Hilbert bases in the hyperplane spanned 
by Hilb(Cio) but this is not true for C 15 . 

Table[5]lists all the 11 tight non-(UHC) cones of dimension 6 that have been found by 
December 17, 2006, including those mentioned already. In the last column we indicate 
whether the Hilbert basis is contained in a hyperplane. 

The smallest non-(UHC), but (ICP) cone we have found has a Hilbert basis of 1 1 el- 
ements, and the largest has a Hilbert basis of 24 elements (most likely (ICP)). Like all 
the others they are extensions of Cio which seems to be the core obstruction to (ICP) and 
(UHC). 

While we have the implications (UHC) (ICP) M = R + M n Z d (under the 
condition that gp(M) = Z d ), it is now clear that the converse implications do not hold. 
However, it remains an open problem whether all cones in dimensions 4 and 5 have (ICP) 
or even (UHC). 
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#Hilb 


#Supp 


ICP 


Aut 


flat 


(-10 


1 n 
10 


2 / 


no 


^20 


yes 




1 o 




no 




yes 


c ! 


1 o 

12 


/I A 

4(J 


yes 




yes 




1 A 

14 




yes 


{id} 


yes 


r' 


14 


39 


no 


{id} 


yes 


r" 


14 


42 


yes 


Z 2 


no 


Q 5 


15 


36 


yes 


{id} 


no 


Cl5' 


15 


36 


yes 


{id} 


yes 




15 


44 


no 


{id} 


yes 


Cl6 


16 


49 


no 


z 2 


yes 


c 


16 


36 


no 


z 2 


yes 



Table 5. Tight non-(UHC) cones 



5. Computational issues 

All programs have been written in C. The tight cones in Section |4] were found by the 
Opteron (O) system mentioned above. It runs Linux, and the executables have been pro- 
duced by the gcc compiler. In the following we will also mention computations on two 
other systems, the author's Intel Core2 6600 (C2) with Windows XP and the DJGPP port 
of gcc, and the University of Osnabriick's Itanium (I) system with Linux and the Intel 
compiler ice. The machines (O) and (C2) are close to each other in speed; (I) is somewhat 
slower (for integer arithmetic), but has very large memory (32 GB). 

Some of the equipment used in the 1998 computations is still accessible. This allowed 
us to measure the gain in speed by improved hardware: the factor is > 40. Moreover, a 
better implementation of shrink yields an acceleration by a factor > 3. In other words, 
4 months of the 1998 search take now a single day. 

The number of cones shrunk by shrink per second depends very much on the pa- 
rameters used for their creation. The performance for 6-dimensional cones generated by 
random 0-1-vectors, whose number varies between 6 and 26, is about 1000 per second 
on (O) or (C2). Cones over 5-dimensional parallelotopes of Euclidean volume < 30 are 
shrunk at a rate of 0.6 per second. The output of tight cones is nevertheless comparable. 

The cones in Section|4]are light food for unicqver. For example, the running time for 
C'( A on (C2) is 1.7 seconds. About 2.5 million vectors are created, but the list of vectors 
in memory simultaneously is bounded by 15,000. 

While all the other programs use 32 bit arithmetic, caradec is set to 64 bit. It has no 
problem with all the cones mentioned, as long as they have (ICP), simply because in all 
cases they are /-covered by cones of determinant < 2. The running time for C" 4 on (C2) 
is 6.9 seconds. 

Despite of its sometimes enormous appetite for memory, caradec has also been suc- 
cessful for all the cones of Table \5\ that lack (ICP), with the exception of Ci6 and C[ 6 . It 
failed for these cones though it was allowed 200 million vectors in memory. 
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The (ICP) property of Ci6 and C[ 6 was falsified by the Monte Carlo method with 1 
million test vectors for each of the non-(UHC) subcones produced by unicover (17.4 
seconds on (C2) for Ck,). 

The longest successful run of caradec with a negative result was C" 5 : 1, 990 seconds 
on (I), 2. 1 billion vectors, 1 10 million simultaneously. The Monte Carlo method does the 
job with 1 million vectors for the single non-(UHC) subcone in 2.7 seconds on (C2). See 
Section [3] for further data C" 5 . 
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